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ABSTRACT 

We investigate the accuracy in the recovery of the stellar dynamics of barred galaxies when 
using axisymmetric dynamical models. We do this by trying to recover the mass-to-light ratio 
(M/L) and the anisotropy of realistic galaxy simulations using the Jeans Anisotropic Multi- 
Gaussian Expansion (JAM) method. However, given that the biases we find are mostly due to 
an application of an axisymmetric modeling algorithm to a non-axisymmetric system and in 
particular to inaccuracies in the de-projected mass model, our results are relevant for general 
axisymmetric modelling methods. We run N-body collisionless simulations to build a library 
with various luminosity distribution, constructed to mimic real individual galaxies, with real- 
istic anisotropy. The final result of our evolved library of simulations contains both barred and 
unbarred galaxies. The JAM method assumes an axisymmetric mass distribution, and we adopt 
a spatially constant M/L and anisotropy (3 Z = 1 — crf/cfj distributions. The models are fitted 
to two-dimensional maps of the second velocity moments V lms = y/V 2 + a 2 of the simula- 
tions for various viewing angles (position angle of the bar and inclination of the galaxy). We 
find that the inclination is generally well recovered by the JAM models, for both barred and 
unbarred simulations. For unbarred simulations the M/L is also accurately recovered, with 
negligible median bias and with a maximum one of just A(M/L) < 1.5% when the galaxy is 
not too close to face on. At very low inclinations (i < 30°) the M/ L can be significantly over- 
estimated (9% in our tests, but errors can be larger for very face-on views). This is in agreement 
with previous studies. For barred simulations the Mj L is on average (when PA= 45°) essen- 
tially unbiased, but we measure an over/under estimation of up to A(M/L) = 15% in our 
tests. The sign of the M/L bias depends on the position angle of the bar as expected: over- 
estimation occurs when the bar is closer to end-on, due to the increased stellar motion along 
the line-of-sight, and underestimation otherwise. For unbarred simulations, the JAM models 
are able to recover the mean value of the anisotropy with bias Af3 z < 0.1, within the region 
constrained by the kinematics. However when a bar is present, or for nearly face-on models, 
the recovered anisotropy varies wildly, with biases up to Af3 z ss 0.3. 

Key words: methods: N-body simulations - galaxies: kinematics and dynamics - galaxies: el- 
liptical and lenticular, cD - galaxies: structure 



1 INTRODUCTION 

The determination of the masses (or equivalently mass-to-light ra- 
tios) of gas-poor galaxies has been an important issue since the 
discovery that galaxies are stellar systems like the Milky Way, 
with mass being a strong driver of many of their properties. Dy- 
namical modelling methods of increased sophistication have been 
developed over the past decades, all based on the assumption 
that galaxies can be described as stationary systems. The first 
attempt at measuring dynamical masses of galaxies were base d 
on the spherical virial equations (|Ppveda 1958]; ISpitzerl Il969l) . 
These methods have the disadvantage that, for accurate results, 
they need to assume self-similarity in the galaxy light and mass 
distribution. More accurate methods allow for axisymmetry and 
take the galaxy light distribution into account. The first d etailed 
axisym metric models of real galaxies were based on the Ijeansl 
equations and assumed a distribution func tion that de- 



pends on two (out of three) integrals of motion (e.g Binnev et al. 
Il99d ; Ivan der Marel et al.lll990l : Emsellem et al.lll994bh . but sve- 
cial classes of three-integral models were also used. Axisymmet- 
ric methods wer e developed to allow for a general orbital distri- 
bution, based on Schwarzschild (1979) numerical orbital superpo- 
sition method (e.g. [C retton et al.lll999l; Iva n der Marel et al II 19981 ; 
IGebhardt et"alll2003l : iThomas et alj|2004 ICappellari et al.N2006l) . 
Currently the most general available models assume galaxies can 
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be approximated by a stationary t riaxial shape (e.g. lde Lorenzi et"al] 
l200llvan den Bosch et all2008h . 

The above modelling techniques were developed under the as- 
sumption that gas-poor galaxies can be well described by stationary 
axisymmetric or triaxial spheroi dal systems. However a key initial 
result of the ATLAS 3D survey dCappellari et alj|2011al hereafter 
Paper I) is the fact that nearby gas-poor gala xies are actually dom- 
inated (86 per cent of them) by fast rotators ( Kr ainovic et alj|2oTll ; 
lEmsellem et alj|201 1[ hereafter Paper II and Paper III), often with 
significant dis k components and resem bling spiral galaxies with the 
dust removed JCappellari et al . 201 lb, Paper VII), 30% of which at 
least are barred. The presence of these bars is a difficult problem for 
all modelling methods and therefore motivates the present study. 

Bars are density waves which results in a tumbling poten- 
tial: this figure rotation is often ignored in the popular dynamical 
modelling methods described above. Dynamical models of barred 
galaxies have been constructe d in the past (e.g. |Pfennigerlll984l ; 
Hafner et al. 200ol; IZhaoll 1996b . However, the existence of intrinsic 
degeneracies in the dynamical modelling of bars make the determi- 
nation of mass quite uncertain even for such models. In fact even the 
full amount of information one can obtain today for external galax- 
ies, namely the full line-of-sight velocity distribution (LOSVD) at 
every position on the sky, is not sufficient to uniquely constrain 
the two free parameters (M/L and inclination) of a simple self- 
consistent axisymmetric mo del ( Vallur fet al.l2004 ; Krainovi c et alj 
l2005l ; ICappellari et al.l2006l : lvan den Bosch & van de Vert2009h . A 
barred model requires at least two extra parameters (the Position 
Angle (PA) and pattern speed of the bar) and dramatically increases 
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the complexity of the orbital structure and the associated degener- 
acy of the problem, instead of improving the accuracy of the mass 
estimate: a broad range of parameters space may well fit the data 
equally well. Moreover, assuming a galaxy is barred also i ncreases 
the d egeneracy in the mass deprojection problem (e.g. Gerhard 
1996), which is al ready mathema tically non unique in the simple 
axisymmetric case (Rybicki 1987). The application of sophisticated 
barred models to large samples would be computationally challeng- 
ing, but feasible exploiting the trivial parallelism of the problem. 
However, this brute-force approach does not remove the intrinsic 
degeneracies so it is not expected to increase the accuracy of the 
mass determinations, and for this reason does not seem justified. 

An alternative approach consists of using some a priori infor- 
mation on the galaxy structure and make empirically-motivated re- 
strictive assumption on the models. This is the approach we are us- 
ing in the systematic determination of the masses of the 260 early - 
type galaxies of the ATLAS 3D survey dCappellari et al.l 120121) . 
We are applying the M ulti-Gaussian Expansion (MGE) technique 
temsellem et al.|[l994al) to accurately describe the photometry of 
all galaxies in the survey (Scott et al. in pre p) and use the Je ans 
Anisotropic MGE (JAM) modelling method dCappellarill2008l) to 
measure masses. The JAM method is based on a simple and very 
efficient solution of the Jeans equations which allows for orbital 
anisotropy ft z = 1 — (<t z /(j_r) 2 . This approach provides good de- 
scriptions of t he integral-field kinematics of the fast rotator early- 
type galaxies dCappellarill2008l : IScott et al.ll2009l ; [Cappellari etail 
2012), which constitute the large majority of the ATLAS 3D sample 
(see Paper II and Paper III). A key motivation for o ur use of the JAM 
method is that it was sh own, using 25 real galaxies dCappellari et al.l 
l2006l : ICappellaril2008l) . to agree well within the model uncertainties 
in the mass determination obtained with the more general axisym- 
metric Schwarzschild approach. 

The use of an axisymmetric dynamical modelling method to 
measure the mass of barred galaxies raises the obvious question 
of what errors in the mass determination are introduced by the ap- 
proach. Answering this question is the goal of this paper. Our work 
is in its spirit an ext ension to barred disk galaxies of the work by 
iThomas et alj d2007l) , which explored the biases introduced by the 
use of axisymmetric models, when extracting masses of triaxial and 
prolate simulated spheroidal galaxy remnants. 

The paper is structured as follows: in Section [2] we briefly 
describe the Multi-G aussian Expansion parametrization from 
lEmsellem et alj dl994al) . the semi-isotropic Jeans equations and the 
method used to create the N-body galaxy models. In Section[3] we 
give an overview of the input models used for the JAM modeling 
in our investigation, and in Section|4]we compare the original and 
recovered values of the corresponding dynamical parameters. Sec- 
tion[5]summarizes all results. 



verts to the mass distribution via the assumed constant M/L. From 
galaxy images this formalism allows us to generate realistic initial 
conditions for our N-body simulations using the method explained 
in Sec 12.31 The MGE parametrization is also the first and crucial 
step of the JAM modeling. Thus, a rigorous and robust approach 
is needed when producing the MGE model of a galaxy, as the pre- 
dicted kinematics may significantly depend on the obtained mass 
distribution (see Sec 12. 2t . The method and softwarq^ we adopt in 
our study t o pro duce MGE parametrization is fully described in 
ICappeilarjd2002h . 

Once the best fit has been found, we have a description of 
the galaxy surface brightness distribution given as a sum of two- 
dimensional Gaussians which we can attempt to deproject. The de- 
projection of a galaxy surface brightness distribution is formally 
non-unique for all but edge-on cases, and the degeneracy can be- 
come severe at low inclinations dRvbickilll987h . The MGE method 
provides just one solution for the deprojection, in terms of a sum 
of three-dimensional Gaussians. This method has been intensively 
used and usually provides luminosity distributions consistent with 
observed photometry of existing galaxies, but the MGE method ob- 
viously does not remove the existing intrinsic degeneracy. 

The deprojection of an MGE mod el can be done anal ytically 
once the viewing angles are known (see lMonnet et al . 1992). When 
the system is assumed to be axisymmetric, only one viewing angle, 
the inclination i (i — 90° for an edge-on system), is sufficient to 
retrieve the full three-dimensional luminosity distribution v (if the 
galaxy is not face-on). 

In a coordinate system (x! , y' , z') centered on the galaxy nu- 
cleus with z' pointing toward us and (x' , y') being the plane of the 
sky, the MGE surface brightness can be written as: 
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where N is the number of adopted gaussian components, each hav- 
ing an integrated luminosity Lk, an observed axial ratio ^ q' k 
1, a dispersion a' k along the major axis, and a position angle (PA) 
ij)k measured counter-clockwise from y' to the major-axis of the 
Gaussian, with (x' r , y' r ) being the associated rotated coordinate sys- 
tem. Then, the deprojected MGE luminosity distribution in cylindri- 
cal coordinates can be expressed as: 
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where the fc-th gaussian has the total luminosity Lk, intrinsic axial 
ratio qk and dispersion a k (In the present study a k = a'k) The 
intrinsic axial ratio qk can then be written as 



/2 2 • 
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(3) 



2 MODELING AND N-BODY SIMULATIONS OF 
EARLY-TYPE BARRED GALAXIES 

2.1 Mass modeling 

2.1.1 Multi-Gaussian Expansion method 

We use i n our study the Multi- Gaus sian Expansion m ethod de- 
scribed in lEmsellem et al.1 dl994al) and lCappellaril d2002l) . The tech- 
nique basically consists of decomposing the luminosity into a num- 
ber of concentric two-dimensional (2D) Gaussians. By fitting the 
detailed surface brightness distribution, the MGE formalism pro- 
vides a description of the intrinsic luminosity density, which con- 



where i is the galaxy inclination. 



2. 1.2 MGE modeling of barred and unbarred SO galaxies 

The fitting of the MGE to the photometry follows the procedure 
applied bv lScott et ai] d2009l) (see their F g. 2) to deal with the pres- 
ence of bars. This same approach is being applied to the MGE fits 
of the ATLAS 3D sample (Scott et al. in preparation). The proce- 
dure allows to find the best fit to the photometry, maximizing the 
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Figure 1. Examples of MGE fits to two model galaxies with a projection angle of i = 25°. The left panel shows the MGE fit to simulated galaxy N4179axi, 
an axisymmetric object. On the right we show the fit to N4754bar with PAj, ar = 60°, where the bar was excluded from the fit. 



minimum q' k and minimizing the maximum q' k , while still being 
consistent with the projected galaxy image within the errors. For 
near face-on cases, a small deviation in the observed axis ratios im- 
plies an important change in the flatness (or roundness) of the mass 
distribution. As the previous procedure frequently ends with very 
similar lower and upper limits for q' k , we force a common axis ratio 
for all Gaussians for such cases to keep an acceptable global shape 
for our MGE model. 



When the bar clearly affects the projected photometry of our 
models (i.e., when it is easily detected), the method adopted for 
the MGE parametrization consists in forcing the lower and upper 
limit of the gaussian axis ratios. In this context, a bar can be con- 
sidered as a perturbation of a disk structure. Bars, if fully fitted by 
MGE components, appear as Gaussians elongated along the appar- 
ent long axis of the bar. The presence of a bar tends to significantly 
affect the q' k values of a few Gaussians depending on its position 
angle, strength and length (the position angle of the bar PAb ar be- 
ing measured counter-clockwise from the projected major-axis of 
the galaxy). The resulting q' k could thus make the system look flat- 
ter (or rounder) than the corresponding axisymmetric ca se if the bar 
is seen end-on (resp. side-on). Previous tests made in IScott et alj 
showed that the best fitting MGE parametrization of a barred 
galaxy is usually not the one which allows the best fit to the ob- 
served kinematics (using JAM models). The kinematic fit is signif- 
icantly improved when the Gaussians have constrained axial ratios 
such that the systems is forced to an axisymmetric "bar-less" MGE 
parametrization. As bars often only affect the photometry within a 
restricted radial range, we use the outer disk of each galaxy to con- 
strain the imposed value of the Gaussian flattening. Figure [T] gives 
two examples of the resulting MGE fits for an axisymmetric simu- 
lation and a barred simulation, both including large-scale disks. 



2.2 Jeans Anisotropic MGE Modeling 

The JAM method is a powerful approach to model the stellar kine- 
matics of early-type galaxies, providing a good description of the 
first two stellar velocity moments (V, K ms ) of a stellar system. 
This technique can be used to probe the dynamical structure of 
ETGs and does in principle allow the recovery of the inclination 
and the dynamical mass-to-light M/L ratio. The JAM technique 
allows for a different M/L and anisotropy for each individual MGE 
Gaussian component. However, the measurement of a global mass 
for real galaxies does not seem to require this extra generality, at 
least with in lR e , where good q uality integral-field data are avail- 
able (e.g. lEmsellem et alj|2004l and Paper I). For this reason the 
models we use make the followi ng simple assum ptions (a full de- 
scription of JAM is provided in Cappellari 2008): 

(i) An axisymmetric distribution of the mass. 

(ii) A constant mass-to-light M/L ratio. 

(iii) A constant anisotropy described by the classic anisotropy 
parameter f3 z = 1 — (<t z /o\r) 2 with f3 z > 0. 

When the mass distribution is represented via an MGE 
parametrization (see Sect. 12.11 above), the Jeans equations 
can be easily integra ted along the line-of-sight as shown by 
lEmsellem et aljl 1994al) in the semi-isotropic case (o z = or ^ <j$) 
and by Cappellaril d2008h in the anisotropic generalization (o z 7^ 
or 7^ (J ti). Here we use t he anisotropic formulas (equations 28 
and 38 of Cappellari (2008)) to derive the projected first and sec- 
ond velocity moments (V and K ms ) given a set of input parameters 
(MGE mass model, mass-to-light ratio M/L, anisotropy f3 z ), and 
thus find the best fitting values within a sampled predefined solution 
space (e.g., /3 Z > 0) 

Such a JAM method has been systemically applied to SAURON 
integral-field stellar kinematics of all 260 early-type galaxies of 
the ATLAS 3D sample. In the present study, we rely on mock ob- 
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Model Galaxy/Model Distance (Mpc) N. M f3 zi j3 z t Os/or time(Gyr) bar 

N4179axi NGC4179 16.5 4e6 0(e) 0.106 1.8 1.5 no 

N4570axi NGC4570 17.1 4e6 0.0 0.145 1.0 1.5 no 

N4442bar NGC4442 15.3 4e6 /3(e) 0.344 1.8 1.5 yes 

N4754bar NGC4754 16.1 4e6 /3(e) 0.343 1.8 1.5 yes 



Table 1. Table providing a list of simulations with their labels and their specifications. Distances are set according to Paper I. N* is the number of particles. /3 z i 
gives the anisotropy of the initial conditions as described in Section l2.3.2l fj z f is the global anisotropy computed for the final state of our simulations. a§ /ctr 
gives the second relation for the geometry of the velocity dispersion ellipsoid for the initial conditions. The time given in the eighth column correponds to the 
simulated time of evolution. The last column indicates whether a bar appeared or not in our simulations. 



Model 


Galaxy/Model 


1 


Dist 


Ns, 


Pzi 


£1 


HernOl 


Hernquist 


1.0 


10.0 


2e6 


0.0 


1.0 


Hern02 


Hernquist 


1.0 


10.0 


2e6 


0.2 


1.0 


Hem03* 


Hernquist(flat) 


0.5 


10.0 


2e6 


0.0 


1.0 


Hern()4* 


Hernquist(flat) 


0.5 


10.0 


2e6 


0.2 


1.0 


N4754ini 


NGC4754 


0.4-0.2 


16.1 


4e6 


0.2 


1.0 



resents a fraction of the total mass, so that given a total number 
of particles per component we can determine Nk, the number of 
particles of that fc-th Gaussian. All components are truncated at a 
chosen radius. To set up the position of each particle, we use a stan- 
dard realisation method with a random generator via the cumulative 
function of a (truncated) Gaussian function, scaling each spatial di- 
mension with the corresponding spatial dispersion. 



Table 2. Table providing a list of the tests simulations (static models) with 
their labels and their specifications. Distance (in Mpc) of iV4754im was 
set according to Paper I, while for Hernquist models it was set arbitrary. The 
flattening of the models is given through the axis ratio q of the MGE models 
(for N4754ini axis ratio cover a range from q = 0.4 for inner gaussians 
to q = 0.2 for outer ones). TV* is the number of particles. j3 z j gives the 
anisotropy of the initial conditions as described in Section 12.3.21 / or 
gives the second relation for the geometry of the velocity dispersion ellip- 
soid for the initial conditions. 



servations computed from N-body simulations of galaxies, and we 
chose to build JAM models from artificial maps to mimic the pro- 
cedure used in the course of the ATLAS 3D survey. We also exclude 
the central few arcseconds during the fitting process, which e.g., 
avoids biases due to the effect of the seeing. We also re-bin all 
m aps before fitting by using the Voronoi tessellation as described 
in ICappellari & Copinl (2003): this allows a guaranteed minimum 
signal-to-noise ratio in each bin, and reduces the scatter in the outer 
parts of the kinematic maps. 



2.3 N-body simulations of regular-rotator galaxies 

As the main motivation of our study is to find the influence of a 
bar on the recovery of basic dynamical parameters with the JAM 
method, we chose to use an N-body approach to generate simu- 
lations of barred early-type galaxies: knowing the exact input dy- 
namics for these simulations, we can then compare the key param- 
eters with those determined via the JAM modeling. We also made 
static realisations of a few (Hernquist and one typical axisymmetric 
lenticular) mass models, and thus only used the initial realisation 
of the N-body distribution. These models are detailed in Table [2] 
while details for evolved simulations are summarized in Table [JJ 
The method to build the initial conditions for our simulations (to be 
evolved, or not) is detailed below. 



2.3.1 Particle positions 

Starting from the MGE parametrization of a mass distribution (after 
taking into account the mass-to-light ratio M/L), the initial posi- 
tions of the particles can be computed easily. Each Gaussian rep- 



2.3.2 Dynamical structure 

Given the particle position, we compute the velocit y dispersio n 
components or, o^ and o z sol ving Jeans Equat i ons (|Jean sl l 19221) , 



within the MGE formalism of lEmsellem et al] dl994ah . For this 
work we use the anis otropic gen e raliza tion of the method (equa- 
tions 19-21 and 34 of Cappellarij ( 1200 81) ), which allows one to set 
arbitrary ratios o z /or and o^/or for the axes of the velocity ellip- 
soid, which is assumed to be cylindrically oriented. Values for these 
ratios can be set individually for each gaussian component, but in 
this study all Gaussians share the same geometry of its velocity dis- 
persion ellipsoid. For some of the simulations initial conditions, we 
used the following (so-called ft — e , with e the intrinsic galaxy 
ellipticity) relation to fix or/o z : 



O.i 



(4) 



This a purely empirical relation which seems to des cribe the general 
trend in the anisotropy of real fast rotator galaxies dCappellari et al.l 
120071) . 



2.3.3 Simulation code 

The numeri cal simulations are pe rformed with a particle-mesh N- 
body code dBournaud et alj|2007l) . The density is computed with 
a Cloud-in-Cell interpolation, and an FFT-based Poisson solver is 
used to compute the gravitational potential, with a spatial resolution 
and softening of 48 pc. Particle motions are integrated with a leap- 
frog algorithm and a time-step of 0.1 Myr. The number of particles 
and the time evolution of each model are given in Table [JJ 



3 INPUTS FOR JAM MODELS 
3.1 Unbarred simulations 

We first built f our simple simulations based on the analytic 
iHernauistl ( fl99(i|) mass distribution 



p{r) 



M a 



2ty r (r + a) 3 



(5) 
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Figure 2. Binned map of the projections of the final state of our simulations with an inclination i = 60° and PAt, ar = 60° for barred simulations. From left 
to right : total projected luminosity / ; mean velocity along the line-of-sight VLos I velocity dispersion along line-of-sight oxos I projected second velocity 
moment V rm s = y ^los °LOS ' P r °fil es °f VL.OS (dot line), ctlos (dash-dot line) and V rm s (solid line) along the galaxy major axis. From top to bottom : 
NA179axi ; N4570axi ; N4A42bar ; A r 47546ar.The square box in the Vrms map and the vertical dot lines in the profiles represent the typical field of view 
used in our study. 



where M is the total mass and a is a scale length. These simulations 
allowed us to quickly check the purely numerical accuracy of the 
JAM modeling method and of our starting-conditions generating 
software for spherical isotropic (HernOl), spherical anisotropic 
(HernOZ), flat isotropic (HernOS) and flat anisotropic (Hern04) 
particle realisations. Flat Hernquist models are oblate systems de- 
rived from the Hernquist profile by forcing an axis ratio of 0.5 
for the mass distribution. For all these simulations the inclination 
was recovered within an error of less than 2° (besides the case 
of HernOl, for which the inclination is meaningless), the global 
anisotropy was accurately recovered within an error of ±0.025 and 
the error on M/L never exceeded 1.5% 

We then considered a more realistic numerical test, using 
the MGE parametrization of the SDSS r-band image of the real 
galaxy NGC4754 and a constant anisotropy and we label this 
model A r 4754im. When computing the particle velocities, we 
forced fi z = 0.2 for all Gaussians. In contrast to the above- 
mentioned Hernquist models, N4754ini represents a complex 
multi-component object in terms of its mass distribution and kine- 
matics, and is therefore expected to be more challenging for the 
JAM modeling method. 



We also built two simulations, respectively based on the MGE 
parametrizations of NGC4570 and NGC4179, which were evolved 
via N-body simulations during 1.5Gyr (the face-on and edge-on 
projections of the final state are illustrated in IA11 . These two 
simulations can be considered as fully relaxed, and contrarily to 
N4754ini and the Hernquist models, the resulting /3 Z is mea- 
sured not to be constant with radius. One important difference be- 
tween N4570axi and N4179axi is that, while the initial condi- 
tions for N4570axi were fixed as isotropic (J3 X — 0), the ones for 
N4179axi were those of a dynamically cold disk as described in 
TableQ] No bar formed either for N4179axi or N4570axi. 

Projected velocity and velocity dispersion maps of our two ax- 
isymmetric simulations with variable f3 z are presented in Figure [2] 
The second velocity moment is dominated by the dispersion in the 
central region, and by the velocity in the outer parts. The profile of 
Kms also shows a central depletion. 

Figure [3] present the local anisotropy j3 z as measured in the 
meridional plane and the equatorial plane, computed on a cylindri- 
cal grid with linearly spaced cells in R, in z and in angle <j> where 
(R,(f),z) are the standard cylindrical coordinates, including a central 
cylindrical cell with a radius of R c — O.Olkpc. This highlights the 
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fact the anisotropy is not constant in our axisymmetric simulations 
and that the central parts are more isotropic than the outer regions. 
In Figure [4] we show that the radial anisotropy profile (which is 
simply the azimuthal average of the equatorial computation) gener- 
ally increases outwards for these axisymmetric simulations. 

3.2 Barred simulations 

We also developed bar simulations, A r 44426ar, and A r 47546ar, 
based on the initial MGE parametrizations of NGC4442 and 
NGC4754, respectively: these are the final state of two N-body sim- 
ulations after 1.5 Gyr of evolution (face-on and edge-on projections 
of the final state can be found in lAU . To generate a bar, we force a 
cold dynamical structure in the initial particle realisations by set- 
ting in the initial conditions aj, = ctr/1.8. The radial velocity 
dispersion ctr was set using the /3(e) function described in Eq. [4] 
With these conditions, a bar appears in each of these simulations 
after only a few rotation period, namely between about 25 Myr and 
50 Myr of simulated evolution. As mentioned, we let the galaxy 
evolve for 1.5 Gyr, to make sure that the bar is well settled. 

N4442bar presents the biggest bar (in size) of our simu- 
lations: we estimate a semi-major axis of 3.0 kpc. The size of 
the bar is determined using the radial flattening (qi SO phote = 
1 - tisov hote) and the position angle (PA iaouhote ) of isophotes 
as done in iMichel-Dansac & Wozniakl {200fj). Basically we define 
the end of the bar as the radius where isophotes are nearly round 
(Qisophote > 0.9) combined with an important change in the po- 
sition angle. Using the same method we determine the semi-major 
axis of the bar for N4754bar to be 2.2 kpc. Outside of the bar re- 
gions, our two simulations are characterized by a rotation pattern 
consistent with an axisymmetric disk-like system. 

The velocity fields of our barred simulations are shown in the 
two lower panels of Figure [2] As for the axisymmetric cases, the 
Vrms maps are dominated by velocity dispersion in the central parts 
and by the mean velocity in the outer parts. But in contrast to the 
N4179axi and iV4754im simulations, the barred simulations all 
show a peak in the center of the Vims maps. This apparent differ- 
ence between the barred and unbarred cases has important conse- 
quences, since this will condition the fit of the projected second 
velocity moment via a JAM model. 

Figure[3]and Figure|4]show that the local anisotropy starts from 
a central value of 0.3, reaching a minimum of nearly 0.2 close to 
the end of the bar, then increasing in the outer parts of the model 
up to 0.7. The derivation of f3 z in the equatorial plane confirms the 
presence of a drop in /3 Z in the outer parts of the bar. 



3.3 Mock observations 

As input for the JAM modeling, we simulated observations by pro- 
jecting our simulations first with four different inclinations: i — 25° 
; i = 45° ; i = 60° and i = 87°, (i = 0° corresponding to the face- 
on projection, and i — 90° to the edge-on projection). The choice 
of a near edge-on projection instead of an exact edge-on projection 
was motivated by the goal of checking the accuracy of the incli- 
nation recovery. Indeed the edge-on projection does not allow an 
overestimation of i and then limits the range of possible uncertainty. 
When a bar is present, we also used four different position angles 
for the bar, and this for each value of the inclination: PAbar = 18° ; 
PAbar = 45° ; PAbar = 60° and PA bar = 87°. The position angle 
of the bar is measured counter-clockwise from the galaxy projected 
major-axis to the bar major-axis, so that PAbar = 87° is close to 
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Figure 4. Azimuthally averaged profile of the anisotropy computed in an 
equatorial plane (— 0.5kpc < z < 0.5fcpc) for axisymmetric simula- 
tions (N4l79axi and N4570axi) and barred simulations (7V44426ar and 
N475Abar). The vertical dash-dot lines correspond to the field of view used 
for the JAM models and therefore the region in which we derived the global 
anisotropy. 
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Figure 5. Effect of the anisotropy parameter f} z on the observed Vrms pro- 
file along major axis at an inclination of i = 18° . The difference between all 
profiles is too small to allow a good recovery of j3 z at such low inclination. 



having the bar end-on and PAbar = 18° close to side-on. We also 
simulated the SAURON pixels size of 0'.'8 assuming a distance for 
each model as given in Table Q] typical of ATLAS 3D objects. Note 
that the M / L provided in Table IC 1 I could have been chosen as unity, 
considering that we will only probe here relative AI/L. However, 
we favoured realistic values as to deal with sensible velocity mea- 
surements. 

For real galaxies the photometry is available over a much 
larger field of view than the kinematics. For each projection, the 
MGE parametrization was achieved using a rather wide field of 
view of 401 x 401 pixels corresponding to 320.8 arcsec, and includ- 
ing the full model. Then, the JAM modeling was fit on a 73 x 73 pix- 
els map (58.4 arcsec) of the second velocity moment with the sim- 
ulated galaxy centered on its nucleus, as this roughly corresponds 
to the setup for a single SAURON exposure. 
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Figure 3. Maps of the local values of the anisotropy f3 z in the equatorial plane and in the meridional plane of N4l79axi (upper left), NA570axi (bottom left), 
A r 4442fear (upper right) and NA754bar (bottom right). Colour gradients for the anisotropy go from f3 z = —0.2 (i.e. <r z > ctr) to /3 Z = 0.8 (i.e. u z < ctr). 
The field of view used in our study is represented by the solid square for the equatorial plane and by the vertical dot line for the meridional plane. 



4 RECOVERY OF PARAMETERS 

The recovery of i, /3 Z and M/L with JAM is done on Kms maps 
only. The prediction of V and a require an extra assumption on 
the constancy of the tangential anisotropy of the JAM models. This 
assumption may not be well verified in the simulations, especially 
in barred ones. However, the accuracy in the determination of the 
above parameters only depends on the ability of JAM to repro- 
duce the Kms, and not the V and a fields se parately. Details o f 
the method to calculate V and a can be found in Cappe llari ( 120081) . 



4.1.1 Anisotropy -inclination degeneracy 



As shown in 



4.1 Recovery of f3 z and inclination 

As shown in Figure[4]the anisotropy varies significantly when going 
from the inner to the outer parts of our simulations. 

The JAM modelling method allows for a different anisotropy 
j3 z for every individual MGE Gaussian. Here we limit ourselves to 
the simple case where f 3 z is constant for the whole model . Tests 
using 25 real galaxies dCappellari etail 120061 : ICappel lari 2008) 
have shown that even with constant anisotropy the recovered M/L 
agrees with the one derived with Schwarzschild models, which al- 
low for a general anisotropy distribution. The extra generality of 
JAM is therefore not required in this case. Thus, for comparison, 
we compute a global anisotropy for our simulations by doing a 
luminosity-weighted average of the local anisotropy for all of our 
simulations. In the following we discuss the issue of the inclination- 
anisotropy degeneracy intrinsic to galaxy dynamics and also the 
important influence of the MGE parametrization on the recovery 
of the global anisotropy. The global anisotropy recovered with the 
JAM modeling method is noted (3 Z AM and the one computed di- 
rectly from the simulation f3 z IM . 



Krai n o vie et aD (2005) and with more galaxies in 



Cappellari et al. I ( feOOot) there is an intrinsic degeneracy in the dy 



namical problem between the recovery of the inclination and the 
anisotropy: given the observed photometry, the observed kinematics 
can be reproduced in detail for a wide range of inclinations, by vary- 
ing the orbital make up of the models. This degeneracy persists even 
in the restricted case in whi ch the anisotropy is assumed to be con- 
stant for the whole galaxy (Cappellari 2008). The degeneracy can 
only be broken by making empirically-motivate d assumpt i ons o n 
the anisotropy. This was the approach adopted by |Cappellaril ( l2008l) , 
who showed using a small sample of galaxies that, if (3 Z is assumed 
to be positiv e, as determined using general models on la rge samples 
of galaxies (Cappellari et al. 20071; Thomas et al . 2009), the correct 
inclination can be recovered from the observed integral-field kine- 
matics. 

Here we test the inclination recovery via JAM using N-body 
simulations of both unbarred (i.e. NAYTdaxi and NA57Qaxi) and 
in particular barred galaxies (i.e. A r 44426ar and N4754bar), for 
which the inclination is known. We confirm the fact that, despite 
the inclination-anisotropy degeneracy the inclination can be recov- 
ered in axisymmetric simulations, and we additionally find that the 
inclination can be recovered even in barred simulations. The influ- 
ence of i on the second velocity moment map is important at low 
inclination. A slight change in the inclination implies a significant 
change in Kms (which is more directly related to the change in the 
axis ratio), and then allows a good recovery of the inclination value 
i. At higher inclination, the change in ellipticity is milder and Kms 
is less influenced by small variations. In all cases, the uncertainty 
never exceeds a few degrees, even when the fit can be considered 
difficult due to the presence of non-axisymmetric features such as 
a bar. We also confirm that at low inclinations (i < 40°) there is a 
degeneracy in anisotropy in the Kms fit, essentially preventing any 
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Figure 6. Recovered values of the global anisotropy for non-barred sim- 
ulations as a function of the projection inclination for N4l79axi and 
N4B70axi. The horizontal lines represent the global anisotropy com- 
puted from the particle realisation (3 z lM (blue for N4179axi and red for 
N4570axi). Vertical lines give the inclination of projections. 



information on the anisotropy from being extracted from the data in 
that case. As shown in Figure [5] changes induced by the anisotropy 
parameter on the V lms profiles are not significant for i = 18°. This 
is due to the fact that at low inclination, maps of the second veloc- 
ity moment are dominated by a z . Therefore, to avoid a too strong 
degeneracy, i — 25° is the lowest inclination we will consider in 
this study. For lower values, the degeneracy prevents f3 z from being 
usefully recovered. 

To summarize, the higher the inclination, the larger the ef- 
fect of the global anisotropy on the projected Kms- For the non- 
barred simulations, following eight different projections (2 galax- 
ies projected at 4 inclinations), we find average errors of 2 de- 
grees and maximum error of 4 degrees at an edge-on view, where 
the deprojected axial ratio changes most weakly with inclination. 
For the barred simulations (32 determinations: 2 galaxies projected 
at 4 inclinations and 4 PAbar) the maximum error in the recov- 
ered inclination is 6 degrees. The maximum error happens only for 
PAbar ~ 90°, when the bar is seen side-on and it is incorrectly 
deprojected as a thin disk (see Fig [8}. The error on pi AM is im- 
portant at low inclination and decreases from face-on to edge-on 
projections. 



4.1.2 Influence of the mass deprojection 

While the inclination-anisotropy degeneracy is intrinsic to the 
galaxy dynamics and determines the uncertainties for the recovered 
global anisotropy and the inclination, /3 Z AM is also influenced by 
the mass deprojection. For regular axisymmetric simulations such 
as N4179axi and N4570axi, the MGE parametrization is rela- 
tively robust, and /?f AM is generally very close to the known intrin- 
sic value /3f IM (see Figure[6}. The deviations near face-on view are 
due to the degeneracy in the mass deprojection at that low inclina- 
tion. 

For the barred simulations N4442bar and N4754bar, the 
presence of the bar implies that axisymmetric models cannot re- 
produce the true stellar density distribution. This means that our 
MGE parametrization can only provide an approximation. Fig [7] 
and Fig [8] illustrate the impact of the degeneracy in the mass de- 
projection of barred galaxies. The deprojected surface brightness 



MGE model of a barred galaxy is shown for three inclination of the 
galaxy (i = 25°, i = 60° and i — 87°) and two positions of the 
bar {PA bar = 18° and PA bar = 87°). Close to edge-on (i = 87°, 
bottom panels), the reconstructed deprojected models do a reason- 
able job at fitting the true edge-on surface bightness contours. When 
the bar is close to end-on, the impact of assuming an axisymmetric 
model becomes more visible. At the other extreme end, near face- 
on models (i = 25°, top panels) have deprojected contours which 
significantly depart from the true edge-on ones. This is mostly due 
to the fact that a small change in the fitted axis ratio of the Gaussians 
has a large impact on the intrinsic axis ratio after deprojection : this 
is further illustrated and emphasised in Appendix IB II For interme- 
diate viewing angles (i — 60°, middle panels), the deprojected pho- 
tometry fits reasonably well the true edge-on contours, while again, 
the discrepancy is emphasised in the region of the bar when it is 
initially viewed edge-on. For a real, observed near face-on galaxy, 
it is hard to know how close the MGE fitting process would get 
from the intrinsic axis ratio of the outer disk, as it would depend on 
e.g., the regularity of the disk (for example, its lopsidedness), the 
signal to noise and contamination from the sky. but Fig|7]and Fig[8] 
demonstrate that a significant degeneracy exists when deprojecting 
targets at low inclination (close to face-on). These tests show that, 
as expected, (i) the presence of a bar produces a deprojected density 
flatter (rounder) that the azimuthally averaged one when the bar is 
close to side-on (end-on), (ii) near a face-on inclination the depro- 
jected density can be significantly in error. This obviously impacts 
the associated dynamical modelling, as detailed in the following 
Sections. 

The results presented in Figure [9] are obtained by applying 
the same method used on observations. We can see that the global 
anisotropy is never well recovered. This is explained by the fact 
that the presence of the bar will tend to flatten or round the Gaus- 
sians in the MGE parametrization if it is seen respectively side-on 
(PAb ar = 18°)or end-on (PAb ar = 87°). It is important to note that 
this bias, being due mainly to the mass deprojection (see Fig|7]and 
FigHJ, is not specific to the adopted dynamical modelling method, 
but is expected to affect more general method like Schwarzschild's 
orbit-superposition method. Moreover the anisotropy is formally 
not a well defined quantity in a barred galaxy, as it is expected to 
vary with the azimuthal location on the galaxy disks. 

To further investigate the impact of the ellipticity e on /?i AM , 
we comput ed the position s of our mock observations in the (V/a, 
e) diagram teinnevll2005l) . (V/a) and e are computed in an ellipse 
of area A — ttR^ where R e is the radius of a cylinder enclosing 
half of the galaxy light. In Figure[lO]we present the result for near 
edge-on projections with R going from 0.5 x i? e to 5.0 x R e . We 
only consider the edge-on case in this figure as all our models fol - 
low the inclination law (see Sec. 4.3 of lBinnev & Tremainell 19871) . 
For comparison additional simple models were also constructed 
and projected edge-on. The first one (named TestOl) is based on 
a MGE parametrization with constant e = 1 — q and isotropic 
kinematics (black solid line on the Figure UOt . The second one 
(Test02) is also isotropic but with a non constant ellipticity in the 
mass distribution, e increasing with radii (black dashed line). 

These two test models highlight the fact that when the con- 
sidered area A is increased (V/a) increases as to roughly follow 
the constant anisotropy lines. The axisymmetric model iV4179a:Ei 
also lay on a constant anisotropy line although its dynamical struc- 
ture presents a j3 z gradient. Barred simulations present a different 
behavior depending on PAbar- With broad field of view (FOV) the 
only effect of the bar is seen through (V/a) which decreases when 
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Figure 7. Effect of the inclination on the deprojection of MGE models for i = 25°, i = 60° and i = 87°. In the first and second columns the MGE fitting 
(red lines) of the projected mass distribution (black contours) is shown for PAi, ar = 18° and PAi, ar = 87°. The edge-on deprojections of the two latter 
MGE fitting are represented in the third and fourth columns respectively, superposed to the azimuthally averaged projected density. 



PAbar increases. But when the size of the FOV is of the order of 
the bar size, e is affected by PAbar and then the projections spread 
over a wide range of anisotropy. The intrinsic ellipticity of the MGE 
parametrization plays here an important role for the anisotropy re- 
covery. 

We made a second test to better understand the effect of the 
MGE parametrization on a model of barred object. From the same 
projection we created two different MGE models : the first one is a 
"free" model with gaussian axis ratios left unconstrained ; for the 
second one we forced the maximum and the minimum axis ratio 
as we do know the intrinsic mass distribution. We did not include 
here the case PAbar = 45° for which the MGE parametrization 
is forced to be neither flatter nor rounder than the axisymmetric 
case. As expected, when we forced the axis ratio during the MGE 
parametrization, the global anisotropy /3f AM recovered was found 
to be much closer to /3f IM as shown in FigQT] We can then assume 
that the accuracy in the recovery of the global anisotropy is mainly 
biased by the MGE parametrization of the photometry. 

Unfortunately, we cannot always objectively choose the best 
deprojected model when we apply the JAM modeling method to 
real observations. And for barred galaxies the accuracy on f3 z will 



only be improved if we can really see the bar or have strong evi- 
dence for its presence. 

4.2 Recovery of M/L 

The Mass-to-Light ratio M/L is the third parameter (after the in- 
clination i and the global anisotropy /3 Z ) computed with the JAM 
modeling method. 

When using the initial conditions, we simply compare the re- 
covered M/L with the input ones. For evolved galaxies we used 
instead, as our reference, the M/L computed from the direct ap- 
plication of the virial relation 2K + W — to the simulation par- 
ticles, where K is the total kinetic energy and W is the total po- 
tential energy of our simulations. The relation makes no other as- 
sumption that a steady state, and thus provides the natural bench- 
mark against which to compare stationary dynamical models. In 
general one expects simulations and real galaxies to satisfy the re- 
lation quite accurately, so that the virial M / L V1I will agree with the 
input one Ml Lsim and no distinction needs to be made. However 
iThomas et alj ( 120071) found that M / L vlr can differ from the input 
one at the 5% level, due to non stationarity, and our results agree 
with theirs. As we are not interested on investigating the stationar- 
ity of the model, but only the biases of the modelling method, for 
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Figure 8. Same as Figure[7]wifh a field of view of 5kpc X 5kpc. 
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Figure 9. Recovered value of the global anisotropy for barred simulations as a function of the projection inclination for A r 4442bar (left) and N47b4bar 
(right) for different PAb ar . The horizontal black solid line represent the value of /3^ IM and the vertical lines shows the inclinations of projection. 
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Figure 10. Positions of edge-on projections of all models in the (V/cr, e) 
diagram. The two tests models described in Section l4.1.2l are in black. The 
green line and the thin black solid lines correspond to isotropy and constant 
anisotropy in the diagram with a step of 0.1 from left to right. The cyan line 
correponds to N4l79axi; the yellow one to A r 4570aici ; the blue ones to 
N4754bar and the red ones to N4A42bar. For barred simulations the solid 
line refers to PAt, ar = 18°; the dashed one to PAt, ar = 45° and the dash- 
doted one to PA bar = 87° (we removed the case PA bar = 60° for better 
legibility). The wider FOV is represented by a square and our reference FOV 
in the present study by a triangle. 




20 30 



Inclination [degrees] 



Figure 11. Recovered values of the global anisotropy as a function of the 
projection inclination for two different MGE parametrization of the same 
model. The color solid lines are the values recovered with a model with 
forced axis ratio for the MGE parametrization and the dash-dot lines are 
for a free MGE parametrization (see Section |4.1.2) . Note that the blue solid 
line and the blue dash-dot line are overlaping. The horizontal black line is 
the value /3^ IM computed from the model. 
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Figure 12. Recovered value of M/L for non barred simulations NAl79axi 
and N4570axi as a function of the inclination of projection. In these regular 
axisymmetric simulations we find an error of 1.5% on M/L (excluding the 
near face-on projection). 
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Figure 13. Recovered value of M/L for A r 47546ar as a function of the 
projection inclination for different values of PA bar . 



maximum accuracy we use as reference M / i v ir in all the compar- 
isons which follow. 

To probe the robustness of our method we applied it to the four 
Hernquist particl e realisations created from the analytic formula of 
Hernquist ( 1990). The results both in the recovery of the global 
anisotropy and the M/L were excellent, with an accurate recovery 
of f3 z and errors on M/L of less than 1.5%. Whilst these simula- 
tions are basic and do not reproduce the complexity of a real galaxy, 
this is a reassuring test of our machinery. 

An intermediate case between simple analytic simulations and 
real galaxies, is the 7V4754im model which is just a regular ax- 
isymmetric rotating galaxy with a constant anisotropy. Its mass dis- 
tribution corresponds to a real galaxy but its intrinsic dynamics are 
simple, as the velocity anisotropy is set to be constant throughout 
the galaxy. The global anisotropy is well recovered within 0.025 
and the mass-to-light ratio is recovered with an error of less than 
1.5%. Although unrealistic, this case helps isolate the influence of 
a variable anisotropy on the results of the JAM modeling method. 

We then used the JAM modeling method on N4179axi and 
N4570axi to explore any systematic bias that may be present with- 
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Figure 14. Recovered value of M/L for A r 44426ar as a function of the 
projection inclination for different values of PAb ar - 



out the presence of a bar. Results are shown in Figure [l2]and JAM 
fitting in Fig. ID II and Fig. ID2I Figure [T2l emphasizes that one can 
expect significant overestimations of M/L up to ~ 10%, when the 
galaxy is close to face- on. This important fact is illustrated with an 
analytic test in Fig. 4 of lCappellari et al. I d2006h and with the galaxy 
NGC0524 in Fig. Al there. Here we confirm it with the present sim- 
ulations. In what follows we will focus on the higher inclinations. 
However, one should keep in mind that the M/L of nearly face-on 
galaxies has to be treated with caution. For i = 45° ; i = 60° and 
i — 87° then, we find that in regular axisymmetric cases M/L is 
recovered with a negligible median bias, and a maximum error of 
just 1.5%. 

Figure [13] and Figure [14] respectively present the results of 
JAM modeling of N4754bar and N4442bar. As for axisymmetric 
cases, we only consider here inclination projections with i > 45°. 
For high inclination projections the M/L is recovered within 3% 
for a PAbar = 45°, which represents the average for random ori- 
entations. However the bias in the recovery can reach up to 15% 
in our tests cases. But the main point is that the recovered M/L is 
correlated with the position angle of the bar, PAb ar . The more the 
bar is seen end-on, the larger the overestimation, due to the larger 
velocities of the stars moving along the bar, and towards the line- 
of-sight, with respect to the other directions in the disk plane. The 
reason is that the presence of a bar produces a peak in the V rms 
maps which is not present in the purely axisymmetric case (see Fig- 
ure[2). In order to fit this peak, the JAM model tend to larger M/L 
values. Moreover, the amplitude of the peak increases with PAbar, 
which explains why the value of M/Ljam also increases with the 
position angle of the bar. The exclusion of the central regions helps 
to reduce the bias introduced by the bar in M/Lj am, but in cases 
where the bar dominates the whole field of view, we cannot expect 
to get rid of its influence. 

The case when the bar is seen side-on (PAbar = 18°) for i = 
45° ; i — 60° and i — 87° is a special configuration for the mass 
deprojection. As the flat bar is deprojected as a flattened disk (see 
FigHJ, following the axisymmetric assumption, the JAM model is 
unable to reproduce the global shape of the Vrms map (see Fig |D4 r 
lD6l Fig lDc^lDTOl top left panel). The error on the recovered mass-to- 
light ratio can be up to ~ 10%. In fact the global shape of the V TU1S 
map is dictated by the global anisotropy. The mass-to-light ratio 
essentially adjusts the fit to the global level of the second velocity 
moment. 
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Figure 15. Recovered value of the mass-to-light ratio as a function of the 
projection inclination. The color solid lines are the values recovered with a 
model with forced axis ratio for the MGE parametrization and the dash-dot 
lines are for a free MGE parametrization. 



The position of the bar PAbar = 45° is a useful reference, as it 
represents the average value for random orientations. As previously 
mentioned in this case the MGE parametrization is hardly affected 
by PAbar- In this configuration the error on M/L does not exceed 
3%, although the Kms map is not reproduced. Basically for pro- 
jections with PAbar < 45° M/L is expected to be underestimated 
while for projections with PAbar > 45° it is overestimated. 

Then for PAbar = 60° and PAbar = 87°, the bar produces a 
vertically elongated structure in projection and an artificially round 
bulge when deprojected as an axisymmetric system. However the 
reproduction of the Kms shape is still a hard task for the JAM 
model (see Fig[D4]-[D6] Fig |P8l IdToI right column). A brief in- 
vestigation pertaining to the influence of the MGE parametrization 
on M/Ljam, illustrated by Figure [15] shows that forcing the flat- 
tening of the mass distribution in the JAM model does not really 
affect M/L} am, except for PAb ar = 87°. For this bar position the 
accuracy in M/Ljam is increased, but at the same time, as pre- 
viously mentioned in Section 147X1 the accuracy on the recovered 
global anisotropy is decreased. The prediction of the V and a fields, 
are however signific antly improved, w hen forcing qu, is as already 
shown in figure 3 of IScott et al] d2009h . 

To sum up, when a bar is present the M / L overestimation or 
underestimation with the JAM modelling method are clearly due : 
(i) to the fact that the mass is incorrectly deprojected as too flat or 
too round, when the bar is edge-on or side-on respectively ; and (ii) 
to the fact that the projected second moments are lower or higher 
when the bar is edge-on or side-on respectively. 

To reduce the effect of the bar, one possible solution is to fit 
the JAM model only to regions where the bar has little or no in- 
fluence. We therefore investigated the effect of the size of the field 
of view (FOV) on M/L jam - One could expect that when increas- 
ing the FOV further we could minimize the effect of the bar on the 
fit and thus reach a better accuracy in the recovering operation. To 
test this, we increased the FOV of our mock data and repeated the 
JAM fitting procedure (the MGE parametrization was not affected 
as it is anyway done on a projection of the simulations with a very 
large field of view). Figure [T6l illustrates M/Ljam as a function 
of the size of one side of the FOV normalized by the size of the 
bar for N4754bar. We find that, as expected, the FOV plays a role 
in the recovered values. When the size of the FOV exceeds the bar 
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Figure 16. Effect of the FOV on the accuracy of the recovered M/L of 
A r 47546ar projected with i = 60° . The x — axis correspond to the size of 
the side of the field of view normalized by the size of the major axis of the 
bar determined in Section f3T2l 

size AM/L decreases and seems to tend to a limit value. This is 
expected as we cannot totally get rid of the effect of the bar in the 
Vrms maps. In our study, the typical fitted field of view is quite com- 
parable to the size of the bar itself, meaning that our previous results 
are close to the worst case scenario. The relative size of the bar with 
respect to the size of the field of view is an important ingredient for 
the recovery of the mass-to-light ratio. 

To conclude, when modeling a barred galaxy assuming an ax- 
isymmetric mass distribution, the M/L is on average (at PAbar = 
45°) still well recovered. It can, however, be overestimated, when 
the bar is parallel to the line-of-sight, or overestimated, when the 
bar is orthogonal to the line-of-sight, by up to 15% (in our tests). 
The amplitude of these errors mainly depend on the position angle 
of the bar PAbar, but also on the size of the FOV. Excluding the 
most central parts and increasing the FOV would naturally tend to 
reduce the influence of the bar without removing the error on M/L 
entirely. 



5 CONCLUSION 

This paper focuses on the study of the possible biases in the Mj L 
and anisotropy determination for barred simulations, when using 
axisymmetric dyna mical models. This extends previous studies 
( Thomas et al. 2007) of triaxial and prolate merger remnants to ax- 
isymmetric and barred simulated disk galaxies that better resemble 
observed fast early- type rotators which constitute the large majority 
of the gas-poor population in the nearby Universe (see Paper II and 
Paper III). We do this by generating N-body simulations of objects 
with properties similar to observed galaxies, both with and with- 
out bars. These are projected at various viewing angles and used to 
generate mock observations that closely resemble real data. These 
data are then fed into the JAM modeling machinery as for real data, 
the difference being that the intrinsic values of the free JAM model 
parameters (i, /3 Z , M/L) are known for the simulated data set. 

The errors in the recovered inclination increase with inclina- 
tion due to the fact that, as previously noticed, the models predic- 
tions are sensitive to the intrinsic axial ratio of the MGE models. 
This implies that for nearly face-on inclinations, where the intrinsic 
axial ratio of the models change rapidly, the inclination is formally 



constrained to a fraction of a degree. This formal accuracy is, how- 
ever, compensated by a broader degeneracy in the mass deprojec- 
tion, leading to a small negative bias in the inclination. In practice, 
for the four simulations we constructed (40 different projections in 
total), the errors never exceeded 5°. 

We confirm previous results that the M/L can be recovered 
within a few percent when the simulated galaxies are nearly ax- 
isymmetric, except for nearly face-on view (i < 30°), for which 
the M/L can be significantly overestimated. The global anisotropy 
Pz can be difficult to recover, especially at low inclination (near 
face-on) due to the inclination-anisotropy degeneracy. This degen- 
eracy implies a significant uncertainty on f3 z at low inclination, but 
a smaller error at high inclination. The global anisotropy is pri- 
marily influenced by the flattening (or the roundness) of the MGE 
parametrization of the projected luminosity. In the case of regu- 
lar axisymmetric objects, the main issue is the intrinsic mathemat- 
ical degeneracy of the luminosity deprojection at low inclination, 
which affects any axisymmetric deprojection method, including the 
adopted MGE one. This results in small deviations of the recovered 
global anisotropy from the value computed from the numerical sim- 
ulations: f3 z is well recovered. 

When a bar is present, the mass deprojection becomes the main 
uncertainty in the models. The deprojected axisymmetric model 
will be naturally different from the true non- axisymmetric barred 
distribution and will change as a function of the observed PA. Con- 
sequently the predicted V^ ms of the models, as well as the corre- 
sponding best fitting f3 z will change as a function of the PA and can 
be quite different from the true axially-averaged value. 

The mass-to-light ratio is less sensitive to the MGE 
parametrization than the global anisotropy, but it is biased due to 
the intrinsic dynamics of the system we want to model. We find that 
M/L is mainly influenced by the position of the bar and the size of 
the field of view. The error depends upon the position angle of the 
bar PAbar and can be up to 15%. Including only regions far from 
the bar allows a reduction of the error, but cannot generally avoid it 
completely. 

Our study provides an estimate of the M/L error that can af- 
fect the determination of dynamical M/L via axisymmetric mod- 
els, and in particular using the JAM method. The large variety of 
possible shapes, sizes and orientations of bars in galaxies, each with 
specific dynamics, prevents us from quantifying the exact errors 
made for individual galaxies. One should also keep in mind that 
our study is done on simulations of relatively weak bars. Therefore, 
the error on the estimated M/L when modelling galaxies exhibit- 
ing stronger bars is expected to be larger. The objects studied here 
are still representative ATLAS 3D fast rotators, this study therefore 
providing clear guidelines when applying axisymmetric modelling 
to such large samples. 
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APPENDIX A: SIMULATIONS PROJECTED DENSITY 

We present here the projected density maps for the simulations used 
in the present study at two different scales. 



APPENDIX B: MGE DEPROJECTION 

We illustrate with the following figure the intrinsic degeneracy 
present in the derojection process of inclined galaxies. 



APPENDIX C: JAM RECOVERY SUMMARY 
APPENDIX D: JAM FITTING 
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Figure Bl. Degeneracy in the edge-on deprojected luminosity with a field of view pf 15kpc X 15fcpc (left column) and 5kpc X 5kpc (right column). The 
first three rows represent the edge-on deprojection of MGE models with an apparent axis ratio q' = 0.907, q' = 0.910 and q' = 0.920 (from top to bottom 
respectively) when i = 25°, supeiposed to the edge-on averaged projected luminosity of NA754bar (black contours). In the last row the projected luminosity 
with i = 25° is ploted for all the three previous MGE models in addition to the luminosity contours of N4754bar with i = 25° and PA\y ar = 87° . We 
see that in this latter case, even though the projected models are indistinguishable, the deprojected ones can be very different. This is a manifestation of the 
unavoidable intrinsic degeneracy in the mass deprojection of axisymmetric bodies at low inclination. 
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Table CI. Table summarizing the values of i, PAb ar , 0z and M/ L of our mock observations and the values recovered by the JAM modelling method. 
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Figure Dl. Comparison between the simulations projected velocity maps and the best JAM fitting for NAl79axi for the four angles of projection. 
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Figure D2. Comparison between the simulations projected velocity maps and the best JAM fitting for NA570axi for the four angles of projection. 
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Figure D3. Comparison between the simulations projected velocity maps and the best JAM fitting for N44A2bar for i = 25° and the four PAb; 
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Figure D4. Comparison between the simulations projected velocity maps and the best JAM fitting for N44A2bar for i = 45° and the four PAb; 
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Figure D5. Comparison between the simulations projected velocity maps and the best JAM fitting for N44A2bar for i = 60° and the four PAb ar - 
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Figure D6. Comparison between the simulations projected velocity maps and the best JAM fitting for N44A2bar for i = 87° and the four PAb; 
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Figure D7. Comparison between the simulations projected velocity maps and the best JAM fitting for A r 4754bar for i = 25° and the four PAb ar - 
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Figure D8. Comparison between the simulations projected velocity maps and the best JAM fitting for A r 4754bar for i = 45° and the four PAb; 
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Figure D9. Comparison between the simulations projected velocity maps and the best JAM fitting for JV47546ar for i = 60° and the four PAb ar - 
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Figure D10. Comparison between the simulations projected velocity maps and the best JAM fitting for A r 4754bar for i = 87° and the four PAb; 
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